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Abstract 

We present here temperature-dependent Raman, x-ray diffraction and specific heat studies be- 
tween room temperature and 12 K on single crystals of spin-ice pyrochlore compound Dy2Ti20j 
and its non-magnetic analogue Lu2Ti20'j. Raman data show a "new" band not predicted by fac- 
tor group analysis of Raman-active modes for the pyrochlore structure in Dy2Ti20'j, appearing 
below a temperature of Tc =110 K with a concomitant contraction of the cubic unit cell volume as 
determined from the powder x-ray diffraction analysis. Low temperature Raman experiments on 
0^^-isotope substituted Dy2Ti20'j confirm the phonon origin of the "new" mode. These findings, 
absent in Lu2Ti20j, suggest that the room temperature cubic lattice of the pyrochlore Dy2Ti207 
undergoes a "subtle" structural transformation near Tc. We find anomalous red-shift of some of 
the phonon modes in both the Dy2Ti20j and the Lu2Ti20i as the temperature decreases, which 
is attributed to strong phonon-phonon anharmonic interactions. 
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I. INTRODUCTION 



Geometrically frustrated magnetic systems are those where spins cannot simultaneously 
minimize all their pairwise exchange energies because of geometrical reasons and hence can- 
not find a unique ground state In a real system, however, any small perturbation to 
these interactions can favor a particular ground state thus lifting the macroscopic degener- 
acy. The perturbation can be in the form of single-ion anisotropy, next or higher neighbor 
interactions, thermal or quantum fluctuations, dipolar interactions, lattice deformations, or 
externally applied stimuli such as pressure or magnetic field. The pyrochlore structure (space 
group Fd3m), with stoichiometry A2B2O7, is a typical example of geometrically frustrated 
magnets in three dimensions. In this structure, A^~^ and B^^ ions form distinct interpen- 
etrating networks of corner-sharing tetrahedra, displaced from each other by half the unit 
cell dimension. If either A or B ion has a magnetic moment, strong geometrical frustra- 
tion results, which in the presence of one or more perturbative terms, gives rise to various 
interesting ground states. Such is the case for the family of insulating rare-earth titanates 
{R2Ti20j] R = trivalent rare-earth ion), which show a wide diversity in their magnetic be- 
havior (For a recent review on these pyrochlores, see ref ^) from a Neel ordered state in 



Er2Ti207 y 



and 



to a spin-liquid state in Tb2Ti207 [51] and spin-ice states in Dy2Ti20j 
Ho2Ti20j |6|. In the spin-ice ground state, two spins on an elementary tetrahedron point 
in and the remaining two point out, an arrangement that is analogous to protons ordering 
in ordinary water ice [3]. Consequently, the ground state degeneracy and hence the "zero- 
point" entropy of a sp.n-ice .s shown B to be very close to what Linus Pauling estimated 
for the water ice nearly 50 years ago [7| . Since spins can easily be manipulated by external 



stimuli, such as magnetic field [9|, |10[ , studying spin ice is of considerable importance even 
from the point of view of understanding the properties of ordinary water ice. For these rea- 
sons, the spin ices Dy2Ti20-j and Ho2Ti20'i have attracted considerable recent attention. 
Besides these two canonical spin ices, other examples of spin-ice or spin-ice like ground state 
include Ho^ Sno Or and Tb^SmOv jl2| members in the analogous stannate series, and. 



Nd2Mo20-j 



13| and Sm2Mo20j ij] in the molybdate pyrochlore series. 



Spin-ice state in Dy2Ti20'j and Ho2Ti20j emerges at very low temperature when the 
dipolar interactions between the large Ising-like Dy^^{Ho^^) spins overwhelm the near- 
neighbor antiferromagnetic exchange, and therefore, a majority of studies on these py- 
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rochlores, until recently, were limited to the liquid Helium or lower temperature range. Re- 
cently. Vibrational properties of some ftanatepyrochlores (includ.ng the spta-i^ Dy.TW.) 



have been studied using infrared reflectivity 



and Raman spectroscopy [16|, |17| up to room 



temperature. In the work by Maczka et al. [l7|] on R2Ti20j (R = Er, Gd and Dy), Raman 
studies of Gd2Ti20j and Er2Ti20-i' were carried out to determine the temperature depen- 
dence of Raman mode frequencies and line-widths. These authors, however, could not follow 
the temperature dependence of Raman modes in Dy2Ti207. These studies have shown that 
mode frequencies of a few phonons decrease upon cooling, a behavior not often observed in 
solids and hence termed as "anomalous" . This anomalous behavior has been qualitatively 
attributed by Bi et al. [15] to strong spin-phonon coupling in Dy2Ti20'j for a few modes. 
Since the strength of exchange interaction is relatively much weaker for titanate pyrochlores 
(ranging between ~ 19 K in TfeaTizOy Q, [l9| to ~ 0.25 K in ^maTzaOy H), and there 
is no long-range spin order in Dy2Ti20-j as well as in Gd2Ti20f and Er2Ti20i down to 
very low temperatures js], this qualitative explanation of the anomalous temperature de- 



pendence in terms of spin-phonon coupling has been questioned |17|. Another explanation 
of this anomalous behavior in terms of a coupling between phonons and crystal field exci- 
tations can also be ruled out because Gd2Ti20-j also exhibit comparable phonon anomalies 
even though in Gd^^ ion no crystal field effect is expected due to spherical symmetry of the 



4f charge cloud. It is suggested in ref. r7| that the anomalous temperature dependence 
could arise from anharmonic interactions. In order to understand the anomalous tempera- 
ture dependence of phonon modes, we have done detailed Raman studies on single crystals 
of Dy2Ti20-j and its non-magnetic analogue Lu2Ti20r over a temperature range of 12 to 300 
K. A quantitative comparison of the phonons of these two pyrochlores is done to understand 
the role of spin-phonon coupling in giving rise to the anomalous phonon behaviors. Most 
interestingly, we show here that an equally anomalous phonon behavior is present for the 
non-magnetic Lu2Ti20j as well. Another important outcome of our study is the observa- 
tion of a "new" mode near 300 cm~^ in the Raman spectra of Dy2Ti20i recorded below T 
= 110 K. Recently, Lummen et al. [l^ also made a similar observation, who tentatively 
assigned the "new" mode to the first excited crystal field level of Dy'^+. By performing 
Raman spectrscopic experiments on an 0^^-isotope enriched crystal of Dy2Ti207, we show 
here for the first time that the "new" mode is actually a phonon mode not expected from the 
factor group analysis of Raman active modes for the space group FdSm of the pyrochlore 
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structure. We also show that the emergence of "new" mode is corroborated by an abrupt 
change in the temperature dependence of the lattice constant by powder x-ray diffraction 
studies. These findings, absent in Lu2Ti'2.0i, suggest a possible lattice symmetry lowering 
in Dy2Ti20j below 110 K. 

II. EXPERIMENTAL TECHNIQUES 

Stoichiometric amounts of Dy20z/Lu20z (99.99 %) and TiOs (99.99 %) were mixed 
thoroughly and heated at 1200 °C for about 15 h. The resulting mixture was well ground 
and isostatically pressed into rods of about 8 cm long and 5 mm diameter. These rods were 
sintered at 1400 °C in air for about 72 h. This procedure was repeated until the compound 
Dy2Ti207 / Lu2Ti207 was formed, as revealed by powder x-ray diffraction analysis, with no 
traces of any secondary phase. These rods were then subjected to single crystal growth 
by the floating-zone method in an infrared image furnace under flowing oxygen. X-ray 
diffraction, carried out on the powder obtained by crushing part of a single crystalline 
sample and energy dispersive x-ray analysis on a scanning electron microscope indicated a 
pure pyrochlore Dy2Ti20-j / Lu2Ti20i phase. The Laue back- reflection technique is used to 
orient the crystal along the principal crystallographic directions. To prepare the 0^^-isotope 
rich samples of Dy2Ti20j, thin crystal slices of Dy2Ti20j (O^^) were cut from their as-grown 
boule. The thickness was further reduced to less than 500 /im by polishing the flat ends. 
Each crystal shce thus obtained was cut into two halves. While one half was annealed in 
oxygen atmosphere enriched with ^^02 (95 %) isotope at 950 °C for 7 days, the other half was 
subjected to an identical heat treatment but under normal oxygen atmosphere. We found 
that the samples annealed under ^^©2 atmosphere gained nearly 2 % in weight, while the 
weight of samples annealed under normal oxygen atmosphere remained unchanged, which 
indicates that a substantial isotopic exchange took place during the heat treatement. 

Raman spectroscopic measurements on a (111) cut thin single-crystalline slice (0.5 mm 
thick and 3 mm in diameter, polished down to a roughness of almost 10 /im) of Dy2Ti20'j 
and Lu2Ti20j were performed at low temperatures in back-scattering geometry, using 514.5 
nm line of an Ar^ ion laser (Coherent Innova 300) with ~ 15 mW of power falling on the 
sample. Temperature scanning was done using a CTI- Cryogenics Closed Cycle Refrigera- 
tor. Temperature was measured and controlled (with a maximum error of 0.5 K) using a 
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calibrated Si-diode sensor and a CRYO-CON 34 temperature controller, respectively. The 
scattered light was collected by a lens and was analyzed using a computer controlled DILOR 
XY Raman spectrometer having three holographic gratings (1800 groves/mm) coupled to a 
liquid nitrogen-cooled charged coupled device, CCD 3000 (Jobin Yvon-SPEX make) of pixel 
resolution 0.85 cm~^. The instrumental broadening is about 5 cm^^. 

High resolution x-ray diffraction measurements were performed between 10-300 K (with 
temperature accuracy better than 0.5 K) using a highly accurate two-axis diffractometer in 
a Bragg-Brentano geometry (focalization circle of 50 /xm) using the Cu-K^ line (A= 1.39223 
A) of a 18 kW rotating anode. 

Specific heat of Dy2Ti20j and Lu2Ti20-j was measured between 5 K and 300 K using 
a Quantum Design Physical Property Measurement System (PPMS). Thin single crystal 
slices (thickness ~ 500 /xm) weighing about 10 mg each were mirror polished and mounted 
on the sample holder of PPMS. Addenda specific heat of the bare sample holder and a small 
amount of APIZONE N grease (to make thermal contact with the sample) was measured 
each time before loading the single crystals on the sample holder. 



III. RESULTS 



A. Raman spectrum of Dy2Ti20j and Lu2Ti207 

The pyrochlore family belongs to the space group Fd3m{0^) with an Aa^sOgO' stoi- 
chiometry, where A^^ occupies the 16d and B'^'^ occupies the 16c positions and the oxygen 
atoms O and O' occupy the 48/ and 86 sites, respectively. Factor group analysis for this 
family of structures suggests six Raman active modes {Aig + Eg + 4^2^) and seven infrared 
active modes {7Fiu). Fig. [T] (bottom panel) shows the Raman spectrum of Dy2Ti20j at 12 
K between wavenumber 50 and 1000 cm~^ which is fitted with fourteen Lorentzians, labeled 
Ml to M14. The assignment of various modes have been done following the previously 



reported results of ref. 2l|, |22] (see Table-I for details). The six Raman active modes are 
labeled M3 (174 cm'^ Fsg), M6 (312 cm'^, Fa^), M7 (330 cm^^ Eg), M8 (453 cm'^ FaJ, 
M9 (515 cm^^, Aig) and MIO (563 cm~^, F2g). The top panel of Fig. [T] shows the Raman 
spectrum of Lu2Ti20-j at 12 K where all the modes are seen, except the mode M5. The six 
Raman active modes of Lu2Ti20'j appear at 188 cm~^ (N3, -^29), 313 cm~^ (N6, -^23), 336 
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(N7, Eg), 458 cm-i (N8,^g), 520 cm-^ (N9, Fa^) and 609 cm-i (NIO, FsJ. Some 



of the previous studies 



23 



24 



25| labeled the mode M9(N9) as Aig + F2g, and overlooked 



Iq have observed a mode near 



the weak mode near 450 cm^^. We and Lummen et al. 
450 cm~^ in several of the pyrochlore titanates. This mode has sometimes been assigned 
to Ti02 as an impurity phase in the pyrochlore sample We rule out this possibility 

because a strong mode of Ti02, near 610 cm~^, is not seen in our sin gle- crystals. We have 



therefore, assigned the 450 cm ^ mode to phonon, as done in ref. 2l|, 1221 ] . 



We now come to the origin of the additional modes M1(N1)^2(N2), M4(N4), Mll(Nll), 



M12(N12), M13(N13) and M14(N14), also observed in refs. 0, ^ and seen in both mag- 
netic and non-magnetic pyrochlores (Fig. [TD . The modes labelled as Mil to M14 can be 

171]. Since the modes Ml, M2 and M4 are seen 



assigned to second order Raman scattering 



in both Dy2Ti20j (magnetic) and Lu2Ti20j (non-magnetic), these can not be assigned to 
crystal field transitions of Dy^^ in Dy2Ti20i. We, therefore, assign these modes to infra- 
red active or silent modes rendered Raman active due to a possible lowering of the local 

n 

symmetry at some crystallographic sites, as suggested in ref. [2^. We have followed the 
evolution of the different modes with temperature from 12 K to 300 K for both Dy2Ti20^ 
and Lu2Ti20j and the results are discussed below. 



B. Evolution of Raman spectra with temperature 

Raman spectra of Dy2Ti20f have been recorded at several temperatures varying from 
12 K to 300 K. The maximum changes are seen in Raman bands Ml to M7 in the spectral 
range of 50 to 400 cm~^. We, therefore, show in Fig. [2] the Raman spectra at a few 
typical temperatures in this spectral range. The modes Ml, M2, M3 and M4 are weak, 
their intensities are, therefore, rescaled as indicated in Fig. [2] with respect to the mode 
M6. It can be clearly seen that for T > 110 K, the spectra in the 250 to 400 cm~^ range 
can be fitted by two Lorentzians representing the two allowed Raman modes (M6 and M7). 
However, the spectra at T = 110 K and below, indicate a "new" mode (termed as M5) 
marked by an arrow in the panels corresponding to T = 110, 50 and 12 K in Fig. [2l The 
frequencies of the modes Ml, M2, M3, M4, M5 (new mode), M6, M7, M9, MIO and Mil 
are plotted against temperature in Fig. [3l It can be seen that the modes Ml, M3, M4, M9, 
Mil and the new mode M5 (which appears below 110 K) show large anomalous red-shifts 
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upon decreasing temperature. The mode M7 initially shows a blue-shift and then a small 
red-shift below 110 K. The mode Ml shows a small blue-shift as temperature decreases below 
110 K. However, the intense mode M6 and the modes M2 and MIO exhibit blue-shift with 
decreasing temperature, an usual behavior seen in solids. The temperature dependence of 
the M3 and M9 modes show an interesting trend, namely, the rate of softening (|^) increases 
below 110 K. 

Temperature dependence of the some of the modes of non-magnetic Lu2Ti20-j is plotted 
in Fig. m It is clearly seen that phonon modes Nl, N2, N3, N4, N6, N7, and N9 show 
anomalous temperature dependence, similar to the ones seen in Dy2Ti20r phonons. 



C. Appearance of a new mode in Dy2Ti20j 

In Fig. [2], it is clearly evident that a "new" mode (M5) emerges near uj ~ 300 cm^^ below 
Tc = 110 K in Dy2Ti20j whereas no new Raman mode has been seen in Lu2Ti20'i even at 



12 K (Fig._^) in the present study or in other titanate pyrochlores, like, Gd2Ti20-r VA , 
Er2Ti207 [ITj, Ho2Ti20j [16| and Sm2Ti20f [20] studied previously. At room temperature, 
the Raman spectrum of Dy2Ti20^ between 250 and 400 cm~^ is satisfactorily fitted to a 
sum of two Lorentzian functions centered at 308 and 322 cm^^, corresponding to the modes 
F2g (M6) and Eg (M7), respectively (see Table-I). A similar fit worked equally well at lower 
temperatures down to about 115 K, as shown for T = 200 and 150 K spectra in the lower 
panels of Fig. [2l However, for the spectrum recorded at = 110 K, it became neccessary 
to include a third Lorentzian in order to capture an essential feature, which is the emergence 
of a "new" mode on the low-frequency side of the mode M6 (middle panel of Fig. [2]). This 
mode gets further resolved and more pronounced as the sample temperature is lowered below 
Tc, as shown by arrows in the panels corresponding to 110, 50 and 12 K. 

We note that Lummen et al. [l^ and M§czka et al. 17] also see this "new" mode in 
their low-temperature Raman spectra of Dy2Ti20-j. Lummen et al. tentatively assigned 
the "new" mode near 300 cm~^ to the first excited CF level after the crystal field splitting 



estimation due to Rosenkranz et al. 28 . However, the crystal field splittings of the lowest 
J-multiplet calculated by Jana et al. 27| differ considerably. Our results of the specific heat 
measurements on Dy2Ti20j and Lu2Ti207, discussed in the later part of this article, suggest 
that the energy gap between the ground and the first excited CF level is ~ 250 cm~^ which 
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is closer to Rosenkranz's estimation than Jana's calculations. Based on these observations, 
the assignment of the "new" mode as a CF mode (transition from ground state to the first 
excited CF level) by Lummen et al. 16| looks attractive. However, this assignment of 
the "new" mode to be the CF origin is not compatible with our low temperature Raman 
experiments on single crystal of 0^^-isotopic Dy2Ti20-j. Since all the Raman active modes 
in pyrochlores involve 0-ion vibrations only, replacement of O^^ by O^^ in Dy2Ti20j will 
red-shift all the phonon modes, whereas the positions of CF modes, if present, will remain 
unchanged. Our experimental results shown in Fig. [5] show a red-shift of ~ 4 cm~^ of the 
"new" mode when O^^ is replaced by O^^, thus suggesting the "new" mode to be a lattice 
vibration and not a CF transition. This assignment is further strengthened by the fact that 
the full width at half maxima (F) is less in O^^ rich Dy2Ti20j (r=6.7 cm~^) as compared 
to its value (F = 9.1 cm~^) in O^^ rich Dy2Ti20T. This decrease in F is also seen for the M6 
and M7 modes. It has been shown earlier that the phonon linewidth depends on isotopic 

|. This, therefore, strongly suggests that the "new" mode M5 near 
300 cm~^ is of phonon origin. Another possibility is that the "new" mode (M5) becomes 
observable below 110 K when the linewidth of the mode M6 decreases. This explanation will 
not be consistent with our results on Lu2Ti20i: at 300 K the linewidth of the N6 mode in 
Lu2Ti20'j is similar to the corresponding M6 mode and no "new" mode is seen in Lu2Ti20-r 
at low temperatures. The appearance of the "new" mode (M5) is, therefore, attributed 
to a Raman inactive phonon mode becoming Raman active due to lowering of the site 
symmetries arising from local structural deformation below 110 K. Following the indication 
of a possible structural transition, we have performed x-ray diffraction experiments as a 
function of temperature and the results are discussed below. 



mass as F cx MJ^^^^^ 



D. Low temperature x-ray diffraction 

In order to examine if the appearance of the "new" mode in Dy2Ti207 at 110 K is related 
to a possible structural transition, we have recorded high accuracy powder x-ray diffraction 
patterns from room temperature down to 10 K. Inspection of the whole diffraction pattern 
recorded at 10 K did not show additional new peaks nor splitting or broadening of the 
existing Bragg peaks. However, as shown in Fig. El the temperature dependence of the 
lattice parameter "a" extracted from (004) and (333) Bragg reflections, clearly displays a 
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discontinuous change in the lattice parameter near 100 K, the temperature at which the 
new mode (M5) appears in our Raman data. In contrast, the temperature variation of the 
lattice parameter in its non-magnetic analogue Lu2Ti20-j, plotted in Fig. E]for comparison, 
is smooth as expected in the absence of a structural deformation. The solid lines in Fig. E] 
are fit to the data by the relation a{T) = adl + ^77^^3112] ? where is the lattice constant 



at K and b and c are fitting parameters 



r(e'=/T'-l)2j' 

29j. For DyiTiiO-j ao=10.116 A, b=3.152 K, 



c=436.3 K and for Lu2Ti20^ they are 10.024 A, 3.119 K and 513.0 K, respectively. 

A survey of the literature for experiments on Dy2Ti20-j showing anomalous behavior near 
110 K directed us to a report 30| on the Nuclear Gamma Resonance of ^^^Dy^^ in Dy2Ti20j 
by Mossbauer spectroscopy suggesting a change in the nuclear quadrupolar interaction with 
temperature. The total electric field gradient in the Dy^~^ nuclei comes from the 4/ orbitals 
of Dy^~^ and the external charges around, in the lattice. As suggested in ref |[30], there is no 
significant change in the lattice contribution with temperature but a dramatic change was 
seen in the electric field gradient due to 4/ orbitals with temperature that becomes almost 
zero at 110 K. We propose a possibility that when the electric field gradient at Dy^^ ions due 
to their 4/ orbitals becomes negligibly small, the surrounding atoms in the lattice slightly 
readjust and move out of their regular Wyckoff positions thus deforming the local crystal 
symmetry without disturbing the average electric field gradient due to the surrounding ions 
in the lattice. There is a recent study on the temperature dependence of spin-relaxation 
time in Dy2Ti20i [31]. Unfortunately, the data are taken only between 10-90 K and at 
300 K. The data corresponding to 300 K do not agree with the value extrapolated from 
the Arrhenius behavior fitted to the data between 10 K and 90 K. Later, experiments by 
Kitagawa et al. ^] filled up the gap in the data in the range of 90 K to 300 K. In Fig. 6 of 



ref. 



321 ]. we can see a small jump near 100 K in the value of the spin relaxation time. It is. 



therefore, tempting to suggest that the discrepancy seen by Sutter et al. |3l| and Kitagawa 



et al. 



321 ] may be related to the structural deformation we observed at 110 K. 



E. Specific heat 

The specific heat of Dy2Ti20-j (Cdto) and Lu2Ti20j {Clto) single crystals in the tem- 
perature range T = 5 to 300 K is shown in Fig. [Tl^a). Specific heat of both compounds 
shows a smooth monotonic decrease upon cooling below room temperature. For Dy2Ti20T, 
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a tendency towards increase in the specific heat values when cooled below 10 K is due to 
short-range magnetic correlations between the Dy'^+ spins jsl- In particular, no anomalous 
feature is seen in Cdto around T = 100 K, the temperature where the "new" mode first ap- 
pears in Raman spectrum. While the specific heat of non-magnetic Lu2Ti20i arises purely 
due to thermal lattice vibrations, the specific heat of Dy2Ti20-j has an additional Schottky 
contribution arising from thermal (de)population of higher lying CF split levels of Dy^+ 
ions. Therefore, Cdto above 10 K can be expressed as an algebraic sum of contributions 
due to Schottky term (Csch) and thermal lattice vibration term (Ciatt)- Since Dy2Ti20j is 
structurally analogous to Lu2Ti20j, albeit with a slightly smaller molecular weight, Ciatt 
term in the specific heat of Dy2Ti20i can be approximated by the measured specific heat 
of non-magnetic Lu2Ti20-j. However, as shown in the inset of Fig. [7](a), Clto in the range 
T = 10 K to 40 K exceeds Cdto by values up to 1 J/R-mol K (~ 10 % near T = 25 K). 
We note that this difference remains large negative even after rescaling Clto by a factor 
= {MdtoI MltoY (fo^ both n = 1/2 and 3/2), where M^o and y^LTO are molecular 
weights of Dy2Ti20f and Lu2Ti20-j, respectively. These observations suggest that the lattice 
specific heat of Dy2Ti20j is smaller than the specific heat of non-magnetic Lu2Ti20-j, which 
is possibly due to the thermal lattice expansion behavior of these compounds shown in Fig. 
ini While the lattice of Lu2Ti20j expands continuously upon heating above 10 K, Dy2Ti20-r 
exhibits a relatively small lattice expansion at first until it expands steeply between 80 and 
100 K and thereafter in a manner analogous to that of Lu2Ti20-j. What transpire from 
the above discussion is that Lu2Ti20i in not a useful lattice template for Dy2Ti207 at low 
temperatures. However, above T = 100 K the thermal expansion behavior and the degree of 
anharmonicity (vide infra) of the two compounds are in quantitative agreement, therefore, 
the specific heat of Lu2Ti20i is used to approximate the lattice specific heat of Dy2Ti20-j 
in the temperature range 100 to 300 K. A rough estimate of the Schottky contribution in 
the specific heat of Dy2Ti20j above T = 100 K is obtained using: Csch{= Cdto ~ Clto) 
■ Csch obtained using this procedure is shown in Fig. Wih) by closed circles. It exhibits a 
broad maximum centered at T = 175 K with a value of ~ 9 J/Dy-mol K. For comparison, 
we^ calculated the Schottky specific heat using the crystal field levels scheme given in ref. 
estimated by rescaling the CF parameters of Ho2Ti20f obtained by Rosenkraz et al. 



28l | using inelastic neutron scattering data. The calculated Schottky anomaly exhibits a 



peak at roughly the same temperature as the experimentally derived Gschi with a value of 
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about 10 J/R-mol K at the peak temperature of 175 K. On the other hand, the calculated 
Schottky anomaly by Jana et al. [27] (shown in Fig. 2 of ref. j27|) exhibits a huge peak 
of 9.5 J/Dy-mol K at a much lower temperature of T = 65 K. Moreover, at still lower 
temperatures, where C^to is less than Clto and hence Csch is expected to be very small, 
the calculated Csch due to Jana et al. ranges from about 1 J/R-mol K near T = 20 K to 
6.5 J/R-mol K near T = 40 K. Our specific heat data clearly favors the crystal field levels 
scheme of Rosenkranz et al., while rule out the validity of the crystal field scheme due to 
271]. From the position of maximum at 175 K in the experimentally derived Cschi 



Jana et al. 



it is easy to appreciate that the first excited crystal field split level of Dy^+ in Dy2Ti20'i lies 
well above room temperature as proposed by Rosenkranz et al. [28]. In Fig. [7|(b) the open 
circles represents the difference Cdto ~ Clto below T = 100 K which disagree with the 
calculated Csch curve, because at low temperatures, use of Lu2Ti20j as a lattice template 
slightly overestimates the Ciatt of Dy2Ti20j leading to even negative values of the difference 
Cdto — Clto below T = 40 K as discussed above. 

In short, the specific heat investigations of Dy2Ti20r and Lu2Ti20i clearly indicate 
that the position of first excited crystal field split doublet is close to 350 K (240 cm~^). 
These results favor the crystal field scheme for Dy^^ in Dy2Ti20-j derived from the work 
of Rosenkranz et al. on Ho2Ti20'! 28|, however, due to problems associated with correct 



determination of Csch at low temperatures, the validity of crystal field splitting due to 
Rosenkranz et al. could not be verified completely. 



F. Anomalous red-shift of phonons 

For spin-ice Dy2Ti20i, as shown in Fig. [3l a large anomalous softening {red-shift) of the 
modes Ml, M3, M4, M9 and Mil is seen with decreasing temperature. The new mode (M5) 
that appears below 110 K also shows a considerable large amount of red-shift with decreasing 
temperature. Interestingly, the phonon modes of the non-magnetic Lu2Ti20i also show large 
anomalous softening, as temperature decreases. As shown in Fig. HJ the phonons Nl, N2, 
N3, N4, N6, N7 and N9 show red-shifts, similar to the Dy2Ti20-j phonons. Temperature 
dependence of the strong F2g modes (M3, N3, M6, and N6) and the Aig mode (M9 and N9) 
of the two pyrochlores are compared in Fig. [8] which emphasizes that magnitude of phonon 
softening with decreasing temperature in both pyrochlores is comparable. 
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Temperature variation of a phonon mode (i) of frequency uoi can be expressed as, 



uj^iT) = u,{0) + (A^,)^,(T) + (A^,)„„,(T) 

+i^^i)ei-AT) + (A.;.).p-,.(T) (1) 

The first term corresponds to the phonon frequency at T = K. The second term is a contri- 
bution to the phonon frequency due to a change in lattice constant with temperature and is 
known as quasi-harmonic contribution to the frequency change. The third term corresponds 
to the intrinsic anharmonic contribution to phonon frequency that comes from the real part 
of the self-energy of the phonon decaying into two phonons (cubic anharmonicity) or three 
phonons (quartic anharmonicity). The fourth term (Aci;j)gj_p^ is due to coupling of phonons 
with charge carriers within the system which is absent in insulating pyrochlore titanates. 
The last term, (Aui) is the change in phonon frequency due to spin-phonon coupling 
arising from modulation of spin exchange integral by the lattice vibration. 

The lattice contribution (Act;j)g^(T) to the phonon energy accounts for an expansion / 
contraction of the lattice leading to a change in the harmonic force constant without changing 
the phonon population. This change in frequency of a mode i can thus be approximately 
related to the change in volume using the Griineisen parameter 7^ = — ^ = ^fpi 
where Bo is the bulk modulus. Taking T=12 K as the reference temperature (the lowest 
temperature in our measurements), '^'^\^^(^^2K)^^ ~ yj]. In order to account 

for this quasi-harmonic contribution, we have taken the measured values of the Griineisen 
parameters from our high pressure Raman and x-ray experiments on Dy2Ti20j at room 



temperature [33] and the values are 7(M3)=2.2, 7(M6)=1.7, 7(M9)=0.9 and 7(M11)=2.2. 
The volume of the unit cell at different temperatures has been taken from the temperature 
dependence of the lattice parameters shown in Fig. O The quasi-harmonic contribution 
is subtracted from the experimental frequencies to obtain the change in frequency (Alj) 
solely due to intrinsic anharmonic effects ((Ac<jj)^^^(T)) and to the spin-phonon coupling 
((AcJi)^p_p;^(T)). This is shown for the modes M3, M6, M9 and Mil in Fig. [91 

The effect of spin-phonon coupling on the phonon frequency may be accounted for using 



15|, 



34 



35j (Au;j)^p_p^(T) = A < SiSj >, where < SiSj > is the spin correlation function 
and A is the spin-phonon coupling coefficient. In the temperature-dependent infrared studies 
of Dy2Ti20j in ref [l^ the phonon anomalies have been attributed to the spin-phonon 
coupling. Our quantitative measurements of the phonon frequency for various modes in 
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the spin-ice Dy2Ti20r (Fig. [8](a)) and the non-magnetic Lu2Ti20j (Fig. [8]^b)) clearly 
show that the magnitude of phonon softening is comparable in both pyrochlores. Since 
Lu2Ti20j is a non-magnetic pyrochlore, there is no contribution from (Ac(jj)^p_p^(T) which 
thus rules out this mechanism for explaining the phonon anomalies in Dy2Ti20i. Therefore, 
the change in frequency (Acu) (as shown in Fig. [HD is solely due to intrinsic anharmonic 
effects ((A^,),„,(T)). 

The anharmonic effects become considerable when the system temperature increases and 
a phonon decay into two or three phonons due to cubic or quartic anharmonic interactions. 
We recall that in the language of quantum many-body theory of phonons in insulating solids 
without isotopic mass disorder are described as quasi-particles with a complex self-energy: 
S = A + iA' originating from anharmonic phonon-phonon interactions. The complex self- 
energy is specific to each phonon and depends on frequency and temperature. The observed 



frequency uJi(T) and damping parameter 7j exhibit a temperature dependence |36l. 1371. l38l|: 



cu,2 = cu,2(0) + 2cu,(0)A(cu,,T) 



^i(O) 



1 + 



A 



and ujiji = 2LJi{0)A'{u!i,T) 



(2) 
(3) 



where a'j(O) is the temperature- independent ith mode frequency in the harmonic approxi- 
mation. Since the highest temperature in our experiment is 300 K, much smaller than the 
Debye temperature of the titanate pyrochlores (~ 500 K), it is sufficient to consider only 
the contribution of cubic anharmonicity. Considering only the simplest decay channels, as 



was done by Klemens 



39| . phonon of frequency uj{q ~ 0) decays into two phonons uJi{qi) 



and co'2((f2) keeping the phonon energy and momentum conserved, i.e. uji{0) = Ui + UJ2 and 
= gi + g2, then. 



AiiiUi, T) = Ai^[l + 2n(cUi(0)/2)] (4) 

where n is the Bose-Einstein mean occupation number of the iih state. The self-energy 
parameter Ai^ can be positive or negative. Most often, the value of Aj"^ is negative, i.e., 
the observed mode frequency decreases with increasing temperature, which is the "regular" 
behavior. However, in case of phonons termed "anomalous" the sign of A/^ has to be 
positive, which is most likely the case for the titanate pyrochlores. Beside the titanate 
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pyroch^ 
Sn02 



ores, other examples of anomalous phonons include: Big Raman-active phonon in 



36|, Big Raman-active phonons in Pri_s5Ceo.i5Cu04^ (see Fig. 7 of ref. |40|), and 



infrared-active TO3 and LO4 phonons in spinels MgAl204 (see Fig. 3 of ref. 4l|) and, the 
one which is most relevant to the present study, F2g{l), i^2g(2) and Eg Raman-active modes 
of /5-pyrochlore K0s20q 42|. Here, the phonon anomaly has been found not only for the 
F2g{X) mode involving potassium, but also for two oxygen modes F2g{2) and Eg, arising 
from highly anharmonic vibrations of oxygen ions. We have discussed the relelvance of this 
to the present case in the next paragraph. It has been argued 38| that the sign of Aj"^ is 
negative or positive depending on whether the restricted two-phonon density of states D{uj) 
(integrated over all decay channels uOi ^ uoi + 002 and = gl + ^2) has a maximum at a 
frequency higher or lower than the phonon frequency. It will be very interesting to obtain 
D{(jj) for Dy2Ti20r from ab initio lattice dynamical calculations to explain the anomalous 
temperature dependance quantitatively. 

We wish to speculate why the phonon-phonon anharmonic interactions are strong in the 
titanate pyrochlores. The most anomalous phonon is the F2g mode near 200 cm~^, which 
is related to O' sublattice vibration. We recall that the 8b site in the pyrochlore structure 
are occupied by O' ions. However, the 8a sites are vaccant. The presence of unoccupied 8a 
sites can result in large vibrational amplitude of O' ions and hence make F2g mode at 200 
cm~^ most anomalous. In order to test the validity of our assumtion, we did a very rough 
comparision of the displacement parameter < > of oxygen ions in /5-pyrochlore K0s20q 
taken from the single- crystal x-ray analysis of ref. j43| with that in the pyrochlore Th2Ti20j 
for which < v? obtained from neutron diffraction data, is available between 4.5 K and 
600 K 4^. We find that the ratio of square- root of vibrational amplitude of oxygen to the 



bond length (Ia-o (A = K or Tb), for the two compounds, have nearly the same magnitude 
(~ 3%). Further, O ion at 48f sites are also slightly displaced from their ideal position |a to 
(| — x)a, where a is the lattice parameter I45I]. The will make the vibrational amplitude of 
O ion at 48f sites larger and hence the Aig mode at 519 cm~^ will also be more anharmonic 
than usual but less than the 200 cm~^ F2g mode involving O' ion. 
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IV. SUMMARY AND CONCLUSION 



In summary, the Raman light scattering properties of the spin ice Dy2Ti20i and its non- 
magnetic analogue Lu2Ti20j have been studied in detail between T = 12 to 300 K. Our 
observations are as follows: (1) In Dy2Ti20j crystal a "new" mode not predicted by factor 
group analysis of Raman-active modes for the pyrochlore lattice emerges upon cooling the 
sample below T = 110 K. Raman data from 0^^-isotopic Dy2Ti20'j strongly indicate the 
phononic nature of this mode. On the other hand, in Lu2Ti20-j, no "new" mode appears in 
the whole temperature range studied; (2) the slope (^) for the modes Ml and M7, shows 
a distinct change of sign at T = 110 K, (3) a large anomalous softening of some of the 
phonons with decreasing temperature is seen in both Dy2Ti20Y and Lu2Ti20j; (4) at the 
same temperature, the cubic lattice parameter of Dy2Ti20-j changes abruptly by about 0.1 
%, while the lattice parameter of non- magnetic Lu2Ti20-j varies smoothly over the entire 
temperature range of our investigations. 

The appearance of a "new" (phonon) mode in the Raman specta is indicative of symmetry 
lowering in Dy2Ti20'j below T = 110 K. However, the powder x-ray diffraction peaks of 
Dy2Ti20i do not show any spliting or emergence of new Bragg peak at least down to 
about T = 10 K. Thus, within the experimental limits of our powder x-ray diffraction 
experiments, the symmetry lowering in Dy2Ti20'i, if any, could not be detected. It should 
be recalled that all six Raman active phonon modes in the pyrochlore lattice are due to 
vibrations of the oxygen ions alone, rare-earth and titanium ions being located at the centre 
of inversion. Since, oxygen ion has a very low x-ray scattering cross-section, a "subtle" 
change in the anionic sublattice may go unnoticed in x-ray powder diffraction experiments. 
We, therefore, propose that a "subtle" structural deformation occurs in the anionic sublattice 
of Dy2Ti20i near = 110 K, giving rise to a "new" phonon mode in the Raman spectra 
below this temperature. A similar situation has recently been noted by Knee et al. |46| for 
the pyrochlore superconductor Cd2Re20'j which exhibits a non-trivial structural deformation 
involving a small deformation in the RcOq octahedral lattice, at 200 K, giving rise to some 
additional low frequency Raman modes. 

Finally, our study of the non-magnetic analogue Lu2Ti20-j which shows an anomalous 
softening of phonons of comparable magnitude as that in Dy2Ti20j suggests that spin- 
phonon coupling is not responsible for the anomalous temperature dependence of phonons 
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in these and other members of the rare-earth titante pyrochlore series. We rather propose 
that strong phonon-phonon anharmonic interaction can lead to these anomalous phonons. 

To conclude, Raman spectroscopy data, in conjunction with x-ray diffraction studies 
of the spin-ice Dy2Ti20j and of its non-magnetic analogue Lu2Ti20i are reported over a 
largely unexplored temperature window, between 10 to 300 K. Our data reveal a "subtle" 
structural deformation of the Dy2Ti20'j lattice at 110 K. We show that anomalous phonon 
softening is a generic feature of rare-earth titanate pyrochlores which is not limited to the 
magnetic members of this family, thereby ruling out the spin- phonon coupling in these 
pyrochlores which was earlier proposed [l5| to be responsible for the anomalous phonon 
behavior. We propose that this behavior is due to unusually strong phonon anharmonicity 
in these pyrochlores. We believe that our experimental results will motivate further studies 
on the spin ice Dy2Ti20-j in order to understand the nature of the phase transition near 110 
K and its possible consequences on the spin ice models. 
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TABLE I: A list of the experimental frequencies in cm at 12 K. 
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"Origin of the mode has been discussed in the text. 



20 




100 200 300 400 500 600 700 800 900 1000 

Raman shift (cm'^) 



FIG. 1: Raman spectrum of Dy2Ti20j and Lu2Ti20r at 12 K. Open circles represent the experi- 
mental data. Raman intensity in the regions from 50-250 cm~^, 425-475 cm"^, 475-650 cm~^ and 
650-1010 cm~^ arc rcscalcd with respect to the region between 250-425 cm~^ in order to view the 
weak modes clearly. The different modes labeled as Ml(Nl) to M14(N14) are assigned in the text 
(Table-I). 
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FIG. 2: (Color online) Raman spectra of Dy2Ti20j at different temperatures. Solid lines are fit 
to the experimental data (open circles), as discussed in text. A new mode (M5) appears at 110 K 
which has been marked by arrows in panels corresponding to 100, 50 and 12 K. Raman intensity 
in the low frequency regions (50-250 cm~^) have been rescaled. 
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FIG. 3: (Color online) Temperature dependence of the different modes of Dy2Ti20j. Solid lines 
are guide to eye. New mode (M5) appears below 110 K 
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FIG. 4: Temperature dependence of the different modes of Lu2Ti207. 
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FIG. 5: Raman spectra of O^^ (open circle) and O^^ (closed circle) isotopic Dy2Ti20'j at 12 K 
in the range of 260 cm-i to 370 cm . A change in frequencies are seen with a change in 0-ion 
isotopic mass. 
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FIG. 6: (Color online) Variation of the lattice parameter of Dy^Fi^O-j with temperature. An 
abrupt jump of the lattice parameter, indicative of a structural deformation, is seen near 100 K, a 
temperature where "new" mode appears in Raman spectrum. In comparison, the lattice parameter 
of non-magnetic Lu^TiiO-j varies smoothly over the entire range of temperature. Solid lines are 
fit to the data as discussed in text. 
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FIG. 7: (Color online) (a) Temperature dependence of specific heat (C) for the pyrochlores 
Dy2Ti20j (DTO) and Lu2Ti20i (LTO) in the temperature range 5 K to 300 K. Inset: blown-up 
view of the low temperature part of the main panel; (b) closed and open circles are data points 
corresponding to the difference Cdto — Clto- The closed circles represent approximate Schot- 
tky contribution (Csch) hi the specific heat of DTO above 100 K. The solid line is the calculated 
Schottky anomaly due to the first five higher lying crystal field doublets at 385, 503, 547, 679 and 
956 K taken from ref. 321 ] . 
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FIG. 8: Percentage change of frequencies (-^%) of (a) Dy2Ti20'j and (b) Lu2Ti20f phonons "with 
temperature "which thus indicates a similar amount phonon softening in both the pyrochlores. 



40 



30 



E 
u 



20 



c 

(0 



s 

<1 



10 



DTO 

• M3 

A M6 

★ M9 
□ Mil 



* * 



★ * 



★ ★ 



★ ★ 



★ * 



★ * 



n . * * " 
□ if 



A A ^ A 



AAA AAA 



^0 100 150 200 250 300 



Temperature (K) 



FIG. 9: The intrinsic anharmonic contribution (AuJanh) to the change in frequency of Dy2Ti20-j 
phonons estimated from experimental data. 
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